EXAMPLES OF SUBSTITUTION SYSTEMS AND THEIR FACTORS 
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Abstract. The theory of substitution sequences and their higher-dimensional analogues is 
intimately connected with symbolic dynamics. By systematically studying the factors (in 
the sense of dynamical systems theory) of a substitution dynamical system, one can reach a 
better understanding of spectral and topological properties. We illustrate this point of view 
by means of some characteristic examples, including a rather universal substitution in one 
dimension as well as the squiral and the table tilings of the plane. 



1. Introduction 

Closed subshifts over finite alphabets are much studied dynamical systems |36} I33j. Their 
understanding in one dimension is a cornerstone of the theory of (symbolic) dynamical sys- 
tems, while systems in higher dimensions show many new phenomena [45]. It is probably 
fair to say that in the latter case, despite great effort, the open questions still prevail. The 
class of symbolic dynamical systems (over finite alphabets) is special also in the sense that 
they possess finitely many (non-periodic) factors (up to isomorphism), where the term 'factor' 
refers to symbolic dynamical systems that are the image of a homomorphism which commutes 
with the shift action (and not to finite subwords). 

Nevertheless, the factors are rarely used explicitly to unravel (details of) the structure of 
the dynamical system. For certain aspects, however, they carry relevant information, and can 
be employed both for structural insight and for concrete calculations. This is particularly true 
for the spectral theory of symbolic dynamics. The latter comes in two flavors. The traditional 
point of view is via the dynamical spectrum [40J, which analyzes the induced shift action on 
the Hilbert space of square integrable functions on the shift space (relative to an invariant 
probability measure). An alternative approach works with the diffraction spectrum [30, 47, 11 1 
of a typical sequence in the space, which was originally motivated by the physical process 
of kinematic diffraction [22J. Its relevance increased with the discovery of quasicrystals by 
Shechtman et al. [46], particularly for systems with pure point (or Bragg) diffraction. Recent 
evidence |48j suggests that also continuous spectral components are practically relevant. 

It is well understood by now that (and how) the two types of spectra are equivalent in 
the pure point case (35j H7J H7], while it has long been known that the dynamical spectrum 
is generally richer in the presence of continuous components \24\ 113] . Nevertheless, in many 
of the classic examples, the complete dynamical spectrum can be reconstructed from the 
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diffraction spectrum of the system and of some of its factors. This can be made more precise 
for systems of symbolic dynamics as well as for Delone dynamical systems with finite local 
complexity [12]. We will see examples of this type below, selected from the class of substitution 
systems with integer inflation factor; see [27l (25J [26], (34] [28] for background material. 

Yet another, and rather recent, aspect is the topological structure of tiling spaces; see 
[44] and references therein. Here, the substitution or the inflation structure can effectively 
be used to calculate the Cech cohomology, for instance via the methods introduced in [7J. 
This provides concrete results, but their (geometric) interpretation is still difficult. This is 
another instance where the careful inspection of the factors can help to understand the detailed 
structure of the cohomology groups. We will explain this in some detail for a recently analyzed 
example, the so-called squiral tiling of the plane [29] HI] . 

Our approach in this article is example-oriented, wherefore we opted for a somewhat infor- 
mal presentation, in the interest of better readability. We assume the reader to be familiar 
with the basic concepts of (symbolic) dynamics, compare [36, 42 j, and with some results from 
their spectral theory, see [40] for more. We begin with an example in one dimension that is 
deliberately designed to have interesting factors, and a rather rich spectrum as a result of this. 
We continue with some remarks on substitution factors with maximal pure point spectrum, 
before we embark on a more detailed analysis of the squiral tiling [291 [14] and its substitution 
factors. This is followed by a similar (though still somewhat preliminary) analysis of the 
classic table tiling [41] . 

2. A FUN EXAMPLE: THE 'UNIVERSAL' MORPHISM 

Undeniably, the best studied substitution rule (on the binary alphabet {a, b}) is given by 
g F : a H > ab, b *— > a and known as the Fibonacci substitution (or morphism). It defines a 
unique hull Xp C {a, 6} z that is strictly ergodic as a dynamical system under the action of 
the (two-sided) shift; see [6] HO] [39] [15] and references therein. It has pure point dynamical 
as well as pure point diffraction spectrum. In a similar way, one can treat all (repetitive) 
Sturmian sequences, a common feature being that they can all be described as model sets 
(also known as cut and project sets); see [37] as well as [39] for a survey. The equivalence of 
the two types of spectra is well understood for systems of this kind [351 El CEB] • 

The relation between the dynamical and the diffraction spectrum is more complicated 
in the presence of continuous spectral components. This was first pointed out in [24] for 
the example of the Thue-Morse substitution, where the dynamical spectrum is richer. In 
particular, a substantial part of the pure point spectrum does not show up in the diffraction 
measure of the Thue-Morse chain, while it can be recovered from the diffraction of the period 
doubling chain. The hull of the latter is a factor of the Thue-Morse system [44 \ 110] . An 
analogous phenomenon was described in [13J for a system with close-packed dimers on Z, and 
is known to also occur for the Rudin-Shapiro chain [15]. The purpose of this section is to 
illustrate these connections with a single substitution rule that entails a mixed spectrum and 
possesses most of the above examples as factors. 
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Let us consider the alphabet {a, b, c,d,a, b, c, d} and define the primitive substitution 

a ab 
b ^ ad 
c i-> cd 
d i-> cb 

which results in a unique hull, for instance via the two-sided sequence 

w = . . . abcbcdadcdcbabcb\abadabcbabadcdad . . . 

which is a fixed point under g 2 with legal seed b\a. The corresponding hull X is obtained as 
the closure of the two-sided shift orbit of w in the product topology. Before we comment on 
the spectral structure, let us identify some relevant factors. 

Upon identifying letters with their barred copies, one obtains the quaternary Rudin-Shapiro 
substitution [51 130] 

£ RS : a i— > ab , b i— > ad, c i— > cd, d i— >• cb , 

with induced fixed point u; Rg of £ RS , with legal seed b\a. This leads to the classic binary 
Rudin-Shapiro chain via the map ip defined by 93(a) = cp(b) = 1 and ip(c) = <p(d) = 1. For 
complexity results, we refer to [5]. The corresponding hulls fail to be palindromic [HE], which 
has interesting consequences for the spectral theory of associated Schrodinger operators. For 
further results on palindromic systems, we refer to [2] and references therein. 

The binary and the quaternary Rudin-Shapiro hulls are mutually locally derivable (MLD); 
see [9] and references therein for the concept. In the symbolic setting, MLD just means that 
there are sliding block maps with local support [34] in both directions of the derivation. The 
non-trivial direction of this claim follows from the observation that 1111 is the longest 1-free 
subword of the binary image and has the unique preimage dcdc under ip. This determines 
the two cosets of Z modulo 2 that are occupied by the letters a, c and 6, d, respectively. Due 
to repetitivity, the subword 1111 occurs with bounded gaps, wherefore this defines a local 
derivation rule. Consequently, both define dynamical systems with the same spectrum, which 
is known to comprise the pure point part Z[^] and an absolutely continuous component, the 
latter being Lebesgue measure with multiplicity 2; see [10] for details. 

Let us go one step back, and define a sliding block map on Xrs that is induced by 

X(ab) = x(cd) = A , x(ad) = x(c6) = B , 
X(ba) = x{dc) = C , xW = x(bc) = D . 

Its action on Xrs is given by u 1— > u x with u x (i) := x(u(i)ii(« + 1)). This way, one induces 
the primitive substitution 

g T : A H> AC , B^ AD, C ■-> BD , D ^ BC , 

which shows a Toeplitz structure (see below) that is somewhat similar to that of the paper 
folding substitution; compare [31 14"! [T5]. 



(1) 



q : 



ai-> ab 
b h- > ad 
c h-> cd 
d h-> cb 
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Let w T be the fixed point of g\ with legal seed C\A, which is the image of w Rg under 
the sliding block map induced by \. This bi-infinite word leads to a partition of Z into four 
subsets A a = {i G Z | w T (i) = a} for a G {A, B, C, D}. The fixed point property then results 
in the equations 

A A = 4Z , A B = 4Z + 2 , 

A c = (8Z + 1) U (16Z + 11) U (4yl c + 3) , 

A D = (8Z + 5) U (16Z + 3) U (4A D + 3) , 

where A A U A B = 2Z and A c U Tl^, = 2Z + 1 was used to simplify the relations. Due to the 
decoupling, one can calculate the solutions via iteration, which gives 

A c = {-1} U U (2 • 4 n+1 Z + (2 • 4 n - 1)) U |J (4 n+1 Z + (3 • 4 n - 1)) , 

n>0 n>l 

A D =\J(2- 4 n+1 Z + (6 • 4 n - 1)) U |J (4 n+1 Z + (4 n - 1)) . 

n>0 n>l 

Note that the singleton set {—1} has to be added to A c because it is not contained in any 
of the lattice cosets (although it is in the 2-adic closure of either of the four unions). The 
only other fixed point of g\ differs from u; T precisely in w T (— 1). It can thus be described 
by moving the singleton set {—1} from A c to A D . Our explicit coordinatization is a result 
of a coincidence in the sense of Dckking |23j. It implies that £ T defines a dynamical system 
with pure point dynamical spectrum, which is Z[^]. The latter coincides with the pure point 
part of the spectrum of the Rudin-Shapiro system. This also follows from standard results; 
compare [23] and references therein. 

Returning to the original substitution g of Eq. ([I]) and mapping all ordinary letters to 1 and 
all barred ones to 1 = —1, one induces the Thue-Morse (or Pruhet-Thue- Morse) substitution 

g TM : 1 i-4 11 , I i-s- 11 , 

formulated on the alphabet {1,1}; compare [B] for background. Our original fixed point 
w is mapped to u> TM = ^tm^tm) with legal seed 1|1. The dynamical spectrum of the 
corresponding hull Xtm under the shift action comprises a pure point part (namely Z[|]) and 
a singular continuous one, the latter leading to an explicit representation as a Riesz product 

Defining the sliding block map ip on Xtm via ip(w)(i) = —w(i)w(i+l), one induces another 
classic substitution, 

g pd : lt-4 11, 

which is known as the period doubling substitution; compare [6]. The image of u> TM under is 
w pd = ^pdC^pd)' w ith legal seed 1|1. A coordinatization of w pd via a partition of Z can be done 
in a similar way as discussed above for g T ; compare [20]. Via the corresponding diffraction 
measure, this provides an alternative way to show that Z[|] is the dynamical spectrum of the 
associated dynamical system. It exhausts the pure point part of the dynamical spectrum of 
the Thue-Morse system. 
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Having analyzed these (selected) factors of our original substitution g, we can conclude 
that the latter has mixed spectrum with all three spectral types being present. The maximal 
equicontinuous (or Kronecker) factor is the dyadic solenoid §3>> which also emerges via the 
'torus parametrization' of the period doubling chain or the maximal model set factor of the 
Rudin-Shapiro chain. Let us explain this type of connection in a little more detail. 

3. Substitution factors with maximal pure point spectrum 

The Thue-Morse dynamical system has a dynamical spectrum of mixed type, with pure 
point and singular continuous components. As was pointed out in [23], the diffraction spec- 
trum only detects part of it. In fact, it only shows the trivial part of the pure point spectrum, 
Z, while the rest is missing (or hidden). However, the period doubling system is a factor that 
is pure point, with dynamical spectrum Z[^], so that the diffraction of Thue-Morse together 
with that of period doubling covers the entire dynamical spectrum of Thue-Morse. 

In fact, the period doubling system is a factor that lies between Thue-Morse and its Kro- 
necker factor, the latter being the dyadic solenoid Si, [23 S3] - This happens to be the 'torus' 
for the period doubling system in the language of model sets [18], and the mapping from 
period doubling to the solenoid is 1-to-l almost everywhere [20} 118]. 

The system of close-packed dimers on the line [13J provides an example where the dynamical 
spectrum has pure point and absolutely continuous components, as does Rudin-Shapiro. In 
both cases, the dynamical spectrum is richer, but the 'missing' parts again can be recovered 
from the analysis of a single factor. A similar situation shows up in many other examples, 
though not in all. Still, it is worth looking at this relation in some more detail. 

A natural generalization of the Thue-Morse substitution, in the spirit of [32] . is 

(k,£) a h-> a k b e 
Q Z™ : b b k a £ 

with k, £ € N, where k = I = 1 is the classic TM substitution. The substitution matrix reads 
\i fc)> w ith eigenvalues k ± I. This two-parameter family of constant length substitutions 
shares many properties with its classic ancestor; see [TO] for a detailed analysis. In particular, 
there is a sliding block map that works for the entire family, and induces a maximal model set 
factor of substitution type. The latter is a generalization of the period doubling substitution, 
namely 

(k,i) . «^ ub 
^gpd : b^ ua 

with the finite word u = b k ~ 1 ab i ~ 1 . This is a constant length substitution with a coincidence 
in the sense of Dekking [23J, and hence a model set [34] . This factor then has a torus 
parametrization [18J which is an almost everywhere 1-to-l map onto a one-dimensional dyadic 
solenoid, which we denote by Sjj, as before. 

Let us now discuss some two-dimensional examples, which are considerably more complex. 
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4. The squiral tiling and its factors 

The squiral tiling appears in [29, Fig. 10.1.4], where it was constructed as a simple example 
for a tiling of the plane by one prototile (and its mirror image) with infinitely many edges. It 
is obtained from the primitive inflation rule 




which has the (linear) integer inflation factor 3. When one distinguishes the two chiralities, 
it is obvious that the (colored) tiling is equivalent (in the sense of mutual local derivability 
|21|, [9] ) to a coloring of the square lattice that emerges from the block substitution 



(3) 



This block substitution is bijective and of constant length in the terminology of [25] ■ In this 
formulation, it was analyzed in detail in [14] . A larger patch is shown in Figure [TJ 

It follows from standard arguments [251 [26] that the dynamical spectrum of the squiral 
contains Z[^] x Z[-|] as its pure point part. Moreover, the analysis of [13] shows that the 
remaining part of the spectrum is purely singular continuous. The corresponding diffraction 
measure is explicitly known, and can be represented as a two-dimensional Riesz product. 

As the tiling is bijective [25], the block substitution ([3]) commutes with the exchange of the 
two colors. This implies that the hull of the squiral is also invariant under the exchange of 
the colors: with every tiling in the hull, also the tiling with inverted colors is in the hull. This 
is the deeper reason why the projection onto the underlying two-dimensional 3-adic solenoid 
S| is at least 2-to-l, so that the (dynamical or diffraction) spectrum cannot be pure point. 
In particular, the squiral tiling cannot originate from a model set. 

Our next step therefore is to find a factor in which pairs of tilings with inverted colors are 
identified. This is achieved with a sliding block map, which maps 2 x 2-blocks of tiles to some 
new symbols. The squiral tiling admits 14 distinct legal 2 x 2-patterns, as only the two blocks 
with four equal symbols are not allowed. In the following block map (in which we write the 
two colors as 1 and 1 = —1), pairs of blocks with inverted colors are identified, but nothing 
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Figure 1 . A square-shaped patch of the squiral tiling with full ^-symmetry, 
as obtained from the block substitution ([3]). For this figure, the reference point 
of the prototiles is placed in their centers, so that the patch is the central part 
of a fixed point under the corresponding block substitution. 
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It is easy to see that any tiling in the image lifts to exactly two squiral tilings. There is one 
position where we can chose whether the lift has a 1 or a 1 there; all other tiles in the lift are 
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then determined. Our factor map is therefore uniformly 2-to-l on the entire hull of all squiral 
tilings. This is completely analogous to the hull of the (generalized) Thue-Morse chains from 
above, which project uniformly 2-to-l to the hull of the (generalized) period doubling chains, 
under a similar sliding block map |10| . 

The block map induces a primitive substitution on the tilings in the factor space, 
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Obviously, the resulting tilings have a Toeplitz structure; compare [3]. They contain lattice- 
periodic subsets, which implies that they are model sets [34J. As such, they necessarily have 
pure point dynamical and diffraction spectra; compare [35]. In turn, the squiral tilings, whose 
hull is a 2-to-l cover of a pure point factor, cannot have pure point spectrum. In fact, with 
suitable scattering strengths (+1 and —1 on the two symbols), their diffraction spectrum is 
purely singular continuous [14] . 

The tilings generated by the substitution ([5|) form the maximal model set factor of the 
squiral tiling, which we henceforth call -F m ax- This factor, in turn, has a two-dimensional, 
3-adic solenoid S| as factor, onto which it projects 1-to-l almost everywhere via the torus 
parametrization of general model sets |18| . We shall now analyze in more detail the set of 
tilings where this projection fails to be 1-to-l. 

For this purpose, we note that the level-n supertiles of the substitution §5§ have the fol- 
lowing block structure 



(6) 



where X is a square block of dimension (3 n — 1) (which is the same for all seven symbols), 
F and G are rows with entries / and g, while E T and G T are columns with entries e and g, 
respectively. 

We can now arrange two supertiles such that the separating line between them passes near 
the origin. Two blocks X are then separated by a line in which all symbols are the same. 
There are four possibilities, two with a horizontal and two with a vertical separation line, 



a h-> 


G 
X 


a 

G T 


5 


b^ 


F 
X 


b 

G T 




C H> 


F 
X 


c 

E T 


5 


d H> 


G 
X 


d 
E T 






e \-± 


G 
X 


e 

E T 


5 




F 
X 


/" 

G T 


5 


g ^ 


G 
X 


g 

G T 









~X~ 




~X~ 


(7) 


F 


1 


G 




X 




X 



X E T X 



X G T X 



FACTORS OF SUBSTITUTION SYSTEMS 



9 



Taking together all ways in which the boundaries of the supertiles can tend to infinity as 
the supertile order is increased, we obtain four one-dimensional model sets, two of which 
are arranged horizontally, and two vertically. Projected to S|, the translation orbits of the 
horizontal pair project to a single translation orbit of a one-dimensional sub-solenoid of type 
S3, and so do the translation orbits of the vertical pair. 

If a quartet of infinite order supertiles is placed such that the common corner remains near 
the origin, the following configurations are obtained, 



(8) 
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Each of these configurations represents a translation orbit in -F max , and all these orbits are 
projected to a single orbit in S|. The configurations shown in (|8|) form the seeds of the seven 
fixed points of the substitution ([5]). 

This structure of the hull is in line with the Artin-Mazur dynamical zeta function [33] of 
the substitution action on the hull of -F m ax- On the one hand, the zeta function is defined in 
terms of the periodic points in the hull under the substitution, 
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where a m is the number of points in the hull that are invariant under an m-fold substitution. 
Likewise, the exponents c m in the Euler product are the cycle numbers, which follow from 
the a m via 

1 , d 



in 



d\m 



where d runs through the divisors of m and jx denotes the Mobius function of elementary 
number theory; see |36l Ch. 6.4] for a detailed exposition. Note that if the hull consists of 
several components for which the periodic points can be counted separately, the total zeta 
function is obtained as the product of the partial zeta functions. In our case, according to the 
analysis above, F m3iX consists of one copy of S§, two extra copies of one-dimensional solenoids 
§3 (above those points where the projection to §3 is 2-to-l), and four extra fixed points above 
the origin of §3. 

This structure of the hull and the dynamical zeta function can be confirmed by computing 
the dynamical zeta function by a second method, which is a by-product of the computation of 
the Cech cohomology of the hull via the Anderson-Putnam complex [7]. The Cech cohomology 
of the hull of a primitive substitution tiling is obtained as the direct limit of the substitution 
action on the Cech cohomology of a finite CW-complex that approximates the tiling space. Let 



10 



MICHAEL BAAKE, FRANZ GAHLER, AND UWE GRIMM 



Table 1. Topology of the squiral hull and its various substitution factors. 
In the left column, we indicate the set on which the projection to 8§ is not 
1-to-l; labels v or h mean that there are vertical or horizontal sub-solenoids 
of type §3 where the projection is 2-to-l, and an integer denotes the extra 
degeneracy (beyond the Id and 2d solenoids) of the projection to the origin 
of S|. In the second column, the cohomology group H 2 is given. In all cases, 
H 1 = Z[g-] 2 and H° = Z. For the factors of F mSLX , examples for identifications 
of the symbols of F max are given in the last column. Especially for the smaller 
factors, there are usually many choices of identifications that yield equivalent 
factors. 
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A( m ) De th e matrix of the substitution action on the m-th cochain group of the approximant 
complex. The zeta function is then given by Thm. 9.1] as 

n fcodd det(i-z^- fc )) _ n fc od d n,(i-A ( ^ fc) ) 



(10) C(z) 



lLevendet(l-^-*)) ^ JJ. (1 _ , A J 



(rf-fc)x 



where the latter equality holds if all ^4( m ) are diagonalizable, with eigenvalues . Note 
that, instead of the action on the cochain groups, one can also take the action on the coho- 
mology (with rational coefficients), as the extra terms in (|10p cancel between numerator and 
denominator. 

For the cohomology of the squiral tiling, we obtain 

(11) H 2 = Z[±] ©Z[i] 2 ©Z 6 , H 1 = Z[\} 2 , H° = Z, 
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while the factor F max leads to 

(12) H 2 = Z[i]©Z[ifeZ 2 ©Z 2 , H 1 = 1[\\ 2 , H° = Z. 

In particular, the cohomology of -F max contains torsion. Taking into account the eigenvalues 
with which the substitution acts on the cohomology, the zeta functions become 

(13) C (z) = (1 = 32)2 = 1 

1 ; ^ sql ' (1 - z)(l - 9z)(l - 3z) 2 (l - z) 3 (l + zf (1 - z)(l - 9*)(1 - z 2 ) 3 ' 

with fixed point numbers Om = 9 m + 4 + 3- (— l) m , and 

CF max (^) 



(1 - z)(l - 9z)(l - 3z) 2 (l - z) 2 (1 - z)3(l - 9*) 

(1 _ 3z) 2 (l-zV 1 



(1 -z)(l -9-z) V 1 - 3 ^/ (i-^) 4 ' 

with c4,f max ^ = 9 m + 3. In the second line of (|14j) ■ we have already expressed the zeta function 
of 

-f^max as a product of the zeta function for S| (where — (3 m — l) 2 ) the square of 
the zeta function for §3 (with d$ = 3 m — 1; compare sequence A024023 of [38]) and the 
zeta function for four extra fixed points. The corresponding zeta functions follow from |16] , 
Alternatively, since solenoids are inverse limit spaces, the cohomology can also be computed 
with the method from [7]. The zeta function of -Fmax clearly reflects the structure of -F max 
already determined earlier. In contrast, the zeta function of the squiral cannot be interpreted 
so easily. 

There is a large number of further factors between F max and the solenoid S§. We have 
systematically studied those factors which can be obtained by identifying some of the seven 
symbols of the substitution ([5]), while still giving rise to a consistent substitution. Such 
identifications of symbols induce identifications of certain tilings projecting to the same point 
on S§. 

Under the identification of symbols / = g, pairs of tilings that project to the same point 
of the vertical sub-solenoid are identified, so that the splitting of preimages above that sub- 
solenoid is closed. What remains is a split horizontal sub-solenoid, and four extra fixed 
points. The cohomology group H 2 loses the torsion part, and one term Z[|] is replaced by Z. 
Closing the splitting of the horizontal sub-solenoid with the identification e = g is completely 
analogous. Under the identification e = / = g, both Id sub-solenoid splittings are closed. 
Only four extra orbits above the origin remain, and H 2 = Z[|] © Z 4 . From here, there are 
many ways to make further identifications. With each such identification, the multiplicity of 
the projection to the origin of S| is reduced by 1, and so is the exponent of the Z fc term in H 2 . 
An example of a smallest factor above S| is obtained with the identifications a = b = c = d 
and e = f = g, with H 2 = Z[§] © Z. Instead of closing the horizontal split solenoid right 
after the vertical one, it is also possible to begin with one or two of the identifications a = b 
and c = d. In such an operation, the exponent of the Z fc term of H 2 drops by 1. Examples 
of MLD classes of the different factors are listed in Table [TJ 
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5. The table tiling 
Infinite legal patterns generated by the block substitution 



(15) 



h-> 



1 

3 



1 ^ 



2 

1 1 



2 



2 1 
2 3 



3 



3 3 
2 



are MLD to the well-known table tilings with mixed spectrum [3T]. The geometric inflation 
rule is explained in Figure [2j and a larger patch is shown in Figure [3j Frettloh |27j has 
independently shown that the table tilings cannot be described as model sets. Indeed, from 
the substitution (fT5|) . one can see that the four supertiles differ in every place, which is a 
generalisation of the bijectivity property. As a consequence, table tilings that consist of a 
single infinite order supertile occur in groups of four, which are pairwise different at every 
position, and which all project to the same point in the underlying solenoid S|. This means 
that the projection to §2 is 4-to-l almost everywhere. In such a situation, the spectrum must 
indeed be mixed, and the tiling space cannot be MLD to the hull of a model set. 



Figure 2. Geometric inflation rule for the table tiling. 




Figure 3. A patch of the table tiling in the form of a level-6 supertile. 
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What can be said about possible model set factors? In the case of the Thue-Morse and 
squiral tilings, we have seen that there exist maximal model set factors, onto which the pro- 
jection has uniform multiplicity. For the table tiling, such a factor with uniform multiplicity 
cannot exist, which can be seen as follows. 

For such a factor to exist, the multiplicity of the projection to S| must be divisible by four 
everywhere. However, there are 10 different pairs of tiles sharing a vertical edge, and 10 pairs 
sharing a horizontal edge, so that for tilings which consist of a pair of infinite order supertiles, 
sharing a vertical or a horizontal line, the projection to §2 ^ s 10-to-l. As a consequence, any 
projection to a model set factor must identify all 10 tilings in such a group of 10, and thus 
cannot have constant multiplicity. 

This does not exclude the existence of model set factors between the table and the solenoid 
§2, however. Such a factor can still have a projection to the solenoid with a non-trivial 
multiplicity at the corners of infinite order supertiles. 

The table tiling admits 24 legal 2x2-patterns, all of which are seeds for one of the 24 fixed 
points under the square of the substitution. These patterns are the following, 



2 
2 



2 

1 



2 
2 1 



1 
3 1 



1 1 

3 3 



1 3 
3 



2 
1 1 



2 1 
1 3 



2 3 




2 3 




2 3 




3 




3 




3 




3 1 




3 3 


2 




1 




2 1 




2 




1 




2 1 




2 3 




2 



2 




1 3 




2 




3 1 




2 




1 




2 1 




3 3 


2 




3 1 




2 




1 3 




1 1 




3 




2 3 




2 



We have found two sliding block maps which induce a primitive substitution on the resulting 
factor. In the first map, the 16 patches of the first two rows map to and the remaining eight 
patches from the third row to 1. This factor map induces a block substitution 



(16) 



£ 

1 



on the alphabet {0, 1}, with I E {0, 1}. The new substitution rule is primitive and admits 
two fixed points (with legal seed \\\\ or \%\\ and reference point in its center, so that the 
fixed point covers Z 2 ). The substitution obviously has a Toeplitz structure, so that each legal 
tiling it generates forms a model set, with pure point spectrum; compare [341 128] . The two 
fixed points project to the same point in S 2 ; so that there are two translation orbits of tilings 
which project to a single translation orbit in S 2 - 

Another (larger) factor is obtained if the first 16 patches are mapped to 0, the first 4 
patches in the third row to 1, and the remaining 4 patches to 2. The induced substitution on 
the alphabet {0, 1, 2} is 



(17) 



i-+ 



"2 


f 


1 H> 


2 


o~ 




2 


2 


1 




1 


2 h-> 


1 




2 


2 


2 



Again, the substitution is primitive and has a Toeplitz structure, implying that also this 
factor is a model set, with pure point spectrum [Ml 28J. There are three fixed points under 
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the square of the substitution, with seeds [ \ 2 ] {t = 1, 2, 3), which project to the same point 
on §!, so that there are three translation orbits of tilings which project to a single translation 
orbit on the solenoid. Everywhere else, the projection is 1-to-l. 

The interpretation of the zeta function turns out to be more difficult if the projection to 
the solenoid fails to be 1-to-l almost everywhere. For the table tiling, we obtain as Cech 
cohomology 

(is) h 2 = z[i]ez[i] 4 ez 3 ez 2 , h 1 = z[±\ 2 , h° = z, 

while the dynamical zeta fuction reads 

(19) Ctable(^) 



(l-4z)(l-4z 2 ) 2 (l-z) 2 (l-z 2 )' 

The corresponding fixed point counts are a m = 4 n + 3 + (— l) n (l + 2 n+1 ). It is natural to 
compare this to the zeta function of 

(1 -2^) 2 

^ = (l-4z)(l-z) ' 

which follows once again from [16] and codes the fixed point counts = (2 n — l) 2 ; compare 
sequences A000225 and A060867 of |38j. At this stage, however, it is far from obvious how 
the fixed point counts of the table shall be mapped to those of the solenoid, under a map 
which is 4-to-l almost everywhere. 
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